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Abstract 

We have studied superfluid-Mott insulating transition of spin-1 bosons interacting antiferro- 
magnetically in an optical lattice. We have obtained the zero-temperature phase diagram by a 
mean-field approximation and have found that the superfluid phase is to be a polar state as a 
usual trapped spin-1 Bose gas. More interestingly we have found that the Mott-insulating phase 
is strongly stabilized only when the number of atoms per site is even. 
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Superfluid-Insulator (S-I) transition has attracted attention, and has been extensively 
studied in the context of 4 He absorbed in the porous media |1|, 0, granular superconductors 
and Josephson-j unction arrays JO. Recently Greiner et al. || has observed an S-I 
transition of 87 Rb atoms trapped in a three-dimensional optical lattice potential by changing 
the potential depth when the number of atoms par site is an integer. This method is the 
most ideal way to study the S-I transition. There are no lattice imperfections, and we can 
easily change the potential depth in order to study both superfluid and insulator phases in 
a single system. 

On the other hand, recent advances of experimental techniques in an optical trap [[|, [?J 
have achieved condensation of spinor bosons. Recent theoretical studies predict a variety 
of novel phenomena of spinor condensates such as fragmented condensation ||, skyrmion 
excitations ||, 10, 11, 12| and propagation of spin waves [12, 13]. 



Here, a natural question is what is expected if the spinor bosons are trapped in an optical 
lattice. In fact, several unique properties of spinor Bose atoms in an optical lattice have been 



suggested by Demler and Zhou [14]. They have proposed some possible phases including the 



superfluid and insulating phases [H], however, no microscopic calculation has been given 
for their boundaries nor stabilities. In this Letter, we study the S-I transition of spin-1 
bosons with antiferromagnetic interaction in an optical lattice at zero temperature when the 



number of atoms per site is an integer. Using a mean- field approximation [115], [TB|, we show 
the zero-temperature phase diagram where the superfluid phase is a polar state as in the 
case of spinor bosons trapped in a usual harmonic trap. More interestingly, we have found 
that the Mott-insulating phase is strongly stabilized only when the number of atoms per 
site is even. 

Bosons with hyperfine spin F — 1, which include alkali atoms with nuclear spin 7 = 3/2 
such as 23 Na, 39 K and 87 Rb, are represented by the Bose-Hubbard model[[n]] in an optical 
lattice as 



<i,j>,a i,a 

+ 2 U ° a l,*,p a lp a iP a i<x + 2 U * S a L a ! T F «/3 ■ F isa i5 a ip , (1) 

i,ct,0 i,a,/3,-f 

where cii a is the annihilation operator for an atom with hyperfine spin a (= 0, ±1) at site i 



and \i is the chemical potential. t=—f drw*(r) (—h 2 V 2 /2M + Vo(r)) Wj(r) is the hopping 
matrix element between adjacent sites i and j, where Wi(r) is a Wannier function localized 
on the ith lattice site, M is the atomic mass and Vo(r) is an periodic potential which 
characterizes an optical lattice. Uq (U 2 ) is the on-site spin-independent (spin-dependent) 
inter-atom interaction. Up (F = 0,2) is defined by Up = cp f c?r|u>j(r)| 4 , where c = 
(go + 2g 2 )/3, c 2 = (g 2 — fl'o)/3, gF = 47ih 2 ap/M, and ap is an s-wave scattering length for 
two colliding atoms with total spin F. We assume an antiferromagnetic interaction U 2 > 
(a 2 > do). This is the case for 23 Na atoms. 

In a Mott insulating phase with large inter-atom interaction (Uo,U 2 3> t), we obtain an 
effective Hamiltonian within the second-order perturbation for the hopping parameter t as 
H ef f = -iiE^Si ■ S,- ./ 2 V. ;/ . (S, ■ S,) 2 = -fE< y >[(S, + S J f-4][(S 1 + S 3 f-4 + 
2/ a], where J x = 2t 2 /(U + U 2 ), J 2 = f + jj^), and a = J 2 / ' J x = U /(U -2U 2 ). 

If we consider the case with only two-sites, we obtain the spin-singlet (highest-spin) ground 
state if U 2 > (<)0. In addition, it is known that at least in one-dimension, the ground state 
of this effective Hamiltonian is a dimerized state with a finite spin excitation gap if U 2 > 0, 



while the ground state is ferromagnetic state U 2 < [18]. The dimerized ground state 
with a positive U 2 suggests the polar state in the superfluid phase as we will see below. 

To study the S-I transition, we use a mean- field approximation in Refs. [|T3| and We 
start from t = case of Eq. [TJ, where the Hamiltonian is reduced to a diagonal matrix with 
respect to sites. Omitting the site index, the single-site Hamiltonian is 



H Q = -^fi + hj fi(h-l) + hj 2 (S 2 -2n), (2) 

where S = a^F^a^ obeys a usual angular momentum commutation relation Sj] = ie^kSk 
and n = Yl a a a a a- S 2 , S z , and n commute with each other. Therefore, the eigenstates of 
the above Hamiltonian are \S, m;n)(— S <m<S), where S 2 \S,m;n) = S(S + l)\S,m;n), 
S z \S,m;n) = m\S,m;n), and n\S,m;n) = n\S,m;n). The energy of the eigenstate is 
E^(S,n) = -fm + \U n(n - 1) + \U 2 [S(S + 1) - 2n\. Since the orbital wave function is 
symmetric, the spin wave function has to be symmetric. As a result, S = 0, 2, 4, n, when 
the number of atoms n is even and 5* = 1, 3, 5, n, when n is odd ||19| . The state with 
m = S is \S,S;n) oc (a\) s (0^)( n ~ s ^ 2 \vac) where 6^ = a^ 2 — 2a 1 ''"a_ 1 "'" creates a spin- 
singlet pair. We obtain other states with lower magnetic quantum numbers by operating 
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S~ = {S + Y to \S,S;n), where S + = S x + iS y = y/2(a^a + a Q ^a_i). Since we assume an 
antiferromagnetic interaction, the ground state is |0, 0; n) with n/2 singlet pairs if n is even, 
whereas the ground state is |1, m;n)(m = 0, ±1) if n is odd. Comparing En \ and 
we obtain the relation between the number of the atoms per site in the ground state 
and the chemical potential (Fig.l). Note that if Uq < 2U2, the atom number per site is even 
all over the phase diagram. 

Let us consider the case of finite t to study the superfluid transition. Here, we introduce 
the superfluid order parameter ip a = (a ia ) = ^/no( a , where no is the superfluid density 
and ( a is a normalized spinor (*( a = 1. The hopping term is decoupled as a\ a aj a ~ 
(■0a a L + ^aUja) ~ ^Pa^a- As a result, the Hamiltonian is represented by a site-independent 
effective Hamiltonian multiplied by the total number of sites. The effective mean-field 
Hamiltonian [16| is 



H mi = Ho + ztJ^^a + V, 

a 

V= -ztyXtp^aq + ijjaal), (3) 

a 

where z is the number of the nearest-neighbor sites. V corresponds to the transfer between 
bosons localized on a particular site and the superfluid. We assume tf) a and t are small and 
include V by a perturbation theory. 

We first consider the case when the number of atoms in a site is even. Calculating the 
second-order correction, we obtain the ground state energy as 

E n (^j) = E(°\0, n)+A (n, U , U 2 , p) • 

1/ n + 3 



A (n,U ,U 2 ,(j) 



1 + 
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3 \fi — U n 

(4) 



—fl + U (n - 1) - 2U 2 

where E = E/zt, p, = /i/zt, Up = Up/zt and if) = (ipi,ip ,ip^i). Since the order param- 
eter is determined to minimize the ground state energy, the ground state is the insulating 
(superfluid) phase with zero (finite) ip if A > (<)0. By the condition A = 0, we obtain the 
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upper (/i + ) and lower (//_) phase boundaries as 

/Z±= -C/ 2 + i[(2n-l)[7o-l] 



± 

6 



^{9f/ 2 + 6(6f/ 2 -2n-3)f7 
+ [36C7| — 12(2n + 3)£7 2 + 9]} 1/2 . (5) 



By equating /i + and we find the minimum of {To, denoting £/q. The result is 



US = -o[QU 2 - (2n + 3)] + -Vn 2 + 3n. (6) 
o o 

The results obtained by the second order perturbation with V cannot determine the 
symmetry of the superfluid order parameter. In order to determine this, we have to calcu- 



late the fourth-order perturbation energy E^> = ^2, npqi Li{i\V\n) 



(n\V\p) (p\V\q) (q\V\t) 
F (0) p (0) „(0) „(0) „(0) „(0) 



E^ 1 Yin V |T(oko P0|- A l° n g but straightforward calculation gives 



E n ^= J B(n,t/ ,t/ 2 ,/i)^|C 2 -2C 1 C-i| 2 

+C(n,U ,U2,n)(fi -J) 2 , (7) 
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with 



B (n, U , U 2 ,fj) 



n(n + 1) 



AE(°)(l,n- l) 2 A£(°)(0,n-2) 

(n + 2)(n + 3) ' 

Af°)(l,n + l) 2 AP(0,n + 2) 



+ 



45 



n(n 



2) 



A£(°)(l,n- l) 2 A£(°)(2,n-2) 

(n + 3)(n + 5) ' 

A£(°)(l,n+ l) 2 A£(°)(2,n + 2) 
n(n + 3) 



15 



1 



A£(°)(l,n + 1) 

n 2 t 



C (n, C/q, U 2 ,\i) 



2 " 



1 



A£?(°)(2,n) 
n(n — 2) 



15 [ AE(°) (1, n - 1) 2 A£(°) (2, n - 2) 

(n + 3)(n + 5) ' 

AK () )(l,n+l) 2 A£( )(2,n + 2) 



n(n + 3) 



45 LAS(°)(l,n 

i n 2 



1 

+ 9 



A£(°)(l,n + i; 
n 



1 



A£(°)(2,n) 

n + 3 



+ 



A£(°)(l,n-1) A£(°)(l,n + 1 
n n + 3 



_A£(°)(l,n- l) 2 + AE(°)(l,n + l) 2 J' ^ 
where AE^ (S, I) = £ (0) (S, l)-E<® (0, n). The first term on the right-hand side of Eq. lifts 
the degeneracy of the superfluid order parameters in spin space. Because B (n, Uq,U 2 ,Ji) 
is negative with even n, the superfluid phase is a polar (spin-0) state ( = (0,1,0) as in 
the spin-1 Bose condensation in a usual harmonic trap [|T^]. The fourth-order calculation 
is explained physically. The second intermediate states are |0, 0; n ± 2), |2, m;n ± 2), and 
|2, m; n). The processes passing through |0, 0; n±2) (|2, m; n±2)) makes negative (positive) 
contributions to i? fn, Uq, U 2 ,pA which favors (disfavors) the polar state. This is physically 
natural since |0, 0; n±2) (|2, m; n±2)) contains a spin-singlet (triplet) pair. Other processes 
passing through |2, m; n) make nontrivial negative contributions to B (n, Uq, U 2 , p) ■ Adding 
all these contributions, we obtain a negative B (n, Uq, U 2 ,Jj), therefore the superfluid is the 
polar state. 

Now consider the case where the number of atoms per site is odd. Since the non- 
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w=3 



U n ! zt 



FIG. 1: Phase diagram of bose-Hubbard model with spin degrees of freedom when U^Jzt = 1. The 
dark region represents the Mott-insulating phase. The dashed lines indicate the phase boundaries 
at t = 0. 

perturbative ground state has degenerated states with m = and ±1, we have to solve 
the secular equation (1, m; n\V E (a)t^ n \ H V\]-) m '] n ) 

to lift the degeneracy. The energy eigenvalues are obtained as 

E® = [-3(0 + 5) - (a + 7 ) + 7(0 + 6)](ft ■ 

± {[ (a _ 7 )_5 (/5 _ 5 )]2 ( ^t.^2 

+4(3/5 + 7 - 25) (a - 2f3 + 36)n 2 \(* - 2CiC i| 2 } V2 , (9) 



where a,/?, 7 , and 5 are given by a = ^ Amo) ^ n _ iy P = V aSW^i) ' T = f Affl)(o,n + i) ' 



and 5 



n+4 



~ Ag(o)(2n+i) res P ectivel y and AE^(SJ) = E^(S,l) - E (0 \l,n). The ground 
state energy corresponds to the lower sign of Eq. |[ Since (3/3 + 7 — 25) (a — 2/3 + 36) is 
positive when n is odd, the ground state of superffuid phase is a polar state as the same as 
the case of even n. 

Surprisingly, the Mott insulating phase with even n is strongly stabilized against the 



superfluid phase comparing with odd-n case. This is intuitively understood as follows: in 
the case of even number of atoms per site, all atoms are able to form singlet pairs on each 
site, while in the case of odd number atoms, one atom remains to be made pairing. In the 
former case, the boson pairs are strongly localized on a site since the formation of singlet 
pairs prevent the bosons from hopping to the nearest-neighbor sites. Since the hopping is 
essential for the superfluid transition, the Mott insulating phase is stabilized in the former 
case. On the contrary, remaining one atom is free to hop to the nearest-neighbor sites more 
freely in the latter case, thus the superfluid transition occurs easily. In addition, Eq. || 
shows that Uq decreases linearly with U2, which is consistent with the above consideration 
since the formation of singlet pairs is energetically more favorable when U2 is larger. 

This "even-odd conjecture" reminds us one-dimensional antiferromagnetic Heisenberg 
models (Haldane's conjecture) []21| and electronic ladder systems |22| such as Hubbard or 
t — J ladders. Both systems show similar properties; the Haldane (ladder) systems have a 
spin excitation gap and an exponential decay of the spin correlation function with an integer 
spin (an even number of legs), while gapless and power- low decay with a half-integer spin 
(an odd number of legs). These properties are able to be essentially explained in terms of 
tightly bound spin singlets as the present study. However, we have applied this even-odd 
conjecture to the Bose systems or S-I transitions for the first time. 

Finally we note that if we assume a ferromagnetic inter-atom interaction ({7 2 < 0), there 
seems to be no strong even-odd dependence of the phase boundaries since the Mott insulator 
phase is the highest spin state and does not include singlet pairs. The detailed comparison 
between the ferromagnetic case and the present study remains as a future problem. Other 
possibilities such as fragmented condensates or two-particle pairings 2^] should also be 
studied in this system. 

To summarize, we have investigated the S-I transition of spin-1 bosons in an optical 
lattice. The zero-temperature phase diagram has been obtained by a mean-field theory. We 
have determined the order parameter of superfluid phase and showed that superfluid phase 
is a polar state. More interestingly, we have obtained a new kind of even-odd conjecture 
that the Mott insulating phase is strongly stabilized against the superfluid phase only when 
the number of atoms per site is even. 
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